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Abstract 



Results are presented from a search for large, extra spatial dimensions in events with 
either two isolated muons or two isolated electrons. The data are from proton-proton 
interactions at ^/s = 7 TeV collected with the CMS detector at the LHC. The size 
of the data sample corresponds to an integrated luminosity of approximately 2 fb^^. 
The observed dimuon and dielectron mass spectra are found to be consistent with 
standard-model expectations. Depending on the number of extra dimensions, the 
limits at 95% confidence level are in the range 2.5 TeV < Ms < 3.8 TeV, where Mg 
characterizes the scale for the onset of quantum gravity. 



Submitted to Physics Letters B 



*See Appendix A for the list of collaboration members 



1 



Models that extend the structure of space-time predict new phenomena beyond the standard 
model (SM) of particle physics. Additional spatial dimensions, essential for formulating quan- 
tum gravity in the context of string theory, have been proposed as a solution to the SM hierar- 
chy problem [I-SI. In this letter, we present a search for events at large dimuon or dielectron 
invariant mass due to contributions from virtual graviton processes in the Arkani-Hamed, Di- 
mopoulos, Dvali (ADD) model [Bill. 

The ADD model postulates the existence of compactified extra dimensions. Gravity is assumed 
to propagate in the entire higher-dimensional space, while particles of the SM are confined to 
a 3-dimensional slice of the multidimensional space. The resulting fundamental Planck energy 
scale Md in the ADD model can be reduced to significantly lower values than suggested by 
the apparent Planck mass Mpi 1.2 ■ 10^^ GeV deduced for 3 spatial dimensions. Md must be 
of the order of the scale of electroweak symmetry breaking to provide an explanation of the 
hierarchy problem. This scenario predicts phenomenological effects that might be observed in 
proton-proton collisions at the LHC. In this paper, we adopt the assumption |4, 5] that all extra 
dimensions are compactified on a torus of size r. In this case, Md is related to Mpi through 
M{?,+2 = M2j/(87rr"),where n is the number of extra dimensions. 

The graviton in this (3 + n) -dimensional formulation can be equivalently expressed as a set of 
3-dimensional Kaluza-Klein (KK) modes [6| with different graviton masses. The coupling of 
the KK modes to the SM energy-momentum tensor leads to an effective theory with virtual- 
graviton exchange at leading order (LO) in perturbation theory. An ultraviolet (UV) cutoff 
A must be introduced to avoid divergences in the summed contributions from all modes. A 
phenomenological consequence of the small mass separation between adjacent KK modes is an 
enhancement in the expected rate of dilepton events at large invariant masses that appears to 
be non-resonant. Depending on the details of the model, virtual-graviton effects can provide 
the dominant experimental ADD signature at high-energy colliders |j4H5|. 

Several ways of parameterizing the LO differential cross sections are provided in the litera- 
ture, including the Han, Lykken, Zhang (HLZ) [4J and the Giudice, Rattazzi, Wells (GRW) [5] 
conventions. In the GRW convention, the leading-order phenomenology for partonic centre- 
of-mass energies ^ A is described by a single parameter Aj, which does not depend on \/s 
for n > 3. The HLZ convention describes the phenomenology in terms of n and a mass scale 
Ms, where Mg is related to the selected UV cutoff and reflects the scale for the onset of quantum 
gravity. Typically, Mg is expected to be of order Md- The parameter Aj can be related to the 
parameters in the HLZ convention fZl. 

The effective theory breaks down at energy scales at which the underlying theory of quantum 
gravity starts to affect the phenomenology. We assume that the range of validity is character- 
ized by a value a/ s^ax, roughly corresponding to the mass Mmax of the lepton pairs emitted in 
the decay of the graviton. As no clear prediction for i/smax can be made within the ADD model, 
and to take into account the requirement \/smax ^ A, most results in this paper are presented 
both for Mmax = Ms and for a range of different values of Mmax- 

Constraints on virtual-graviton signatures in the ADD model of extra dimensions have been 
obtained at HERA [8, 9], LEP IITOlfell. and the Tevatron dlllZl- At the LHC, limits have been 
presented based on measurements of diphoton events l|T8l - |20| . 

CMS uses a right-handed coordinate system with axes labeled x, y, and z, and the origin at 
the center of the detector. The z-axis points along the direction of the anticlockwise beam. 
The azimuthal and polar angles are (p and 6, with 9 measured from the positive z-axis. The 
pseudorapidity rj is defined hy t] = — lntan(0/2). 
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A main feature of the CMS apparatus is a superconducting solenoid of 6 m internal diameter, 
providing a magnetic field of 3.8 T. Located within the field volume are silicon pixel and strip 
inner trackers, an electromagnetic calorimeter (ECAL), and a hadronic calorimeter (HCAL). 
The ECAL consists of lead-tungstate crystals covering pseudorapidities of |?/| < 1.5 (barrel) 
and 1.5 < \f]\ < 3.0 (endcaps). The CMS muon detectors are embedded in the return yoke 
of the magnet. Muons are measured with detection planes using three different technologies: 
Drift Tubes, Cathode Strip Chambers, and Resistive Plate Chambers. The first stage of the CMS 
trigger system employs custom hardware and processes information from the calorimeters and 
the muon system. The event rate is further reduced by a computer farm using the event infor- 
mation from all detector systems. A detailed description of CMS can be found in Ref. ||2T|. 

This analysis uses data samples collected with the CMS detector in 2011, corresponding to an 
integrated luminosity [22] C of 2.3 ± 0.1 fb^^ (dimuons) or 2.1 ± 0.1 fb^^ (dielectrons). The 
integrated luminosity for the dimuon channel is larger because the muon selection has less 
stringent requirements on the performance of the calorimeters during data-taking. The event 
selection follows closely the criteria used in a search for Z' bosons ||23| . The muon data sample 
was collected using a single-muon trigger with a transverse momentum (pj) threshold which 
was varied between 15 and 40 GeV over the course of data-taking to allow for changes in in- 
stantaneous luminosity. The selection of electron events is based on a trigger requiring the 
presence of 2 electrons or photons with energy depositions > 33 GeV. Candidate events are 
required to have a reconstructed interaction vertex with |z| < 24 cm, and a radial distance 
y^;c^ + 1/2 < 2 cm. For events passing the complete selection requirements, the trigger efficien- 
cies for signal and SM Drell-Yan (DY) events with large mass are > 99%, with an uncertainty 
of < 1%. 

Muons with \ri\ < 2.1 and px^^, > 45 GeV are selected. The candidates are required to be identi- 
fied both in the outer muon system and the inner tracker, and the inner track must contain re- 
constructed energy deposits in at least 1 pixel layer and more than 10 strip-tracker layers. Muon 
candidates are required to have signals from at least two muon detector layers included in the 
reconstructed muon track. Muon candidates satisfy the isolation requirement Pj/ Pxpt < 0.1, 
where the sum extends over the momenta p'j of all charged particle tracks (excluding the muon 

/ 2 2 

track) within a cone of size AR = -J {Atjy + {A(py = 0.3 around the muon direction of flight. 
To reject backgrounds from cosmic-ray muons, we require a transverse impact parameter rel- 
ative to the primary vertex of < 0.2 cm, and an opening angle of (Xpipi < n — 0.02 between the 
2 muon momentum vectors. No charge requirement is applied to the muon pairs. However, 
all selected muon pairs of mass > 450 GeV are found to have opposite charges. Events with 2 
muons passing the selection criteria are accepted for analysis. 

Electron candidates are reconstructed from energy depositions in the ECAL (superclusters) 
matched to a track in the silicon tracker. ECAL superclusters are constructed from 1 or more 
clusters of energy depositions surrounding the crystal with the highest local energy deposition. 
An associated track is required to contain signals from at least 5 tracker layers. The track must 
be matched geometrically to the supercluster, and the spatial distribution of energy must be 
consistent with that expected for an electron. Only electron candidates with a ratio of energy 
depositions in the HCAL and ECAL below 0.05 are considered. To minimize the contamination 
from jets, electron candidates are required to be isolated. Candidates with a sum of transverse 
track momenta > 5 GeV within 0.04 < AR < 0.3 around the candidate track are rejected. In 
the ECAL and inner HCAL layer, the deposited transverse energy in a cone AR = 0.3 
around the electron candidate (excluding the transverse energy Ej^e of the electron) must be 
< 2 GeV + 0.03 x £T,e for the barrel, or < 2.5GeV(<2.5GeV + 0.03 x {E^, - 50 GeV)) for the 
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endcaps if Ej^e < 50GeV(£T,e > 50GeV). Additionally, the Ej deposition in the outer HCAL 
layer within 0.15 < AR < 0.3 around the electron position is restricted to < 0.05 GeV. Selected 
events must contain 2 electrons of opposite charge, each with transverse energy Ej > 35 GeV 
(in the barrel region), or with Ej > 40 GeV (in the endcaps). The explicit charge requirement is 
found to have negligible influence on the presented results. Events in which both electrons are 
reconstructed in the endcaps are not used in the analysis, since electrons from the ADD signal 
would on average be produced at smaller values of tj than the SM backgrounds. 

The search is performed with a set of events that contains either electron or muon pairs above 
a mass value Mmm- The lower bound of the signal region is chosen to maximize the expected 
upper limits of the ADD model parameter Aj in each lepton channel. The optimum value of 
Mmin is found to be 1.1 TeV for both the dimuon and the dielectron channel, based on simulation 
studies. 

In both search channels, the PYTHIA 8.142 HMUFl event generator with the MSTW08 [26] par- 
ton distribution function (PDF) set is used to simulate the expected signal. Interference terms 
between the standard model DY process and the virtual graviton are taken into account in the 
evaluation of the signal cross sections. Simulated events for both signal and SM backgrounds 
are passed through a detailed detector simulation based on GE ant4 |27|, using a realistic CMS 
alignment scenario, and the same reconstruction chain as data. 

In this analysis. The SM DY process is the dominant background. In the dimuon channel, we 
use the MC@NLO ||28l |29| event generator with the CTEQ6.6 L30J PDF set to simulate the DY 
background. The parton level events from MC@NLO are passed to HERWIG 6 ||3l1 for the sim- 
ulation of the QCD parton shower and hadronization, PHOTOS {32] for the simulation of the 
electroweak (EW) parton shower, and Jimmy [33J for the simulation of multiple parton interac- 
tions. The simulated reconstruction efficiencies in the chosen region of acceptance, including 
all selection criteria, are found to be 90% ± 3% for the high-mass DY dimuon background and 
90% ± 4% for the ADD dimuon signal. 

Mass-dependent corrections [34J beyond the QCD next-to-leading order (NLO) predictions im- 
plemented in MC@NLO are studied to improve the SM DY estimate in the dimuon channel. EW 
NLO effects are evaluated by comparing HORACE [35J NLO predictions interfaced to HERWiG 6 
with HORACE LO predictions interfaced to HERWIG 6 and PHOTOS. In this procedure, PHOTOS 
corrections are applied to the LO results to account for radiation effects as these corrections 
are also included in the DY simulation based on MC@NLO. The effect of electroweak NLO 
corrections is found to be smaller than the QCD NLO contribution and of opposite sign. The 
estimated correction factor for the DY background beyond 1.1 TeV is ~ 0.90 ± 0.06. Next-to- 
next-to-leading order (NNLO) QCD corrections are obtained using higher-order calculations 
from FEWZ [36 J. The corresponding multiplicative correction factor for the DY background in 
the signal region is estimated to be 1.03 ± 0.03. For the purpose of setting limits, both the EW 
NLO and QCD NNLO corrections are applied to the DY background prediction obtained from 
MC@NLO. 

The DY background in the dielectron channel is simulated with PYTHIA 6 137| and normalized 
according to the observed data in the range 60-120 GeV around the Z resonance. As in the 
dimuon channel, electroweak NLO corrections at large mass are studied with HORACE. The 
estimated correction factor for the DY background beyond 1.1 TeV is found to be 0.92 ± 0.06. 
The simulated reconstruction efficiency for the high-mass DY dielectron background and the 
ADD dielectron signal in the selected acceptance range, including all selection criteria, is found 
to be 84% ± 3%. 
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The parton distribution functions have a strong impact on the SM DY background in both 
search channels. The PDF uncertainties for the DY process are evaluated by comparing results 
from the CTEQ6.6 [30J, MSTW08 [26\, and NNPDF21 [38J PDF groups. This procedure follows 
the recommendations of the PDF4LHC working group ||39| . The uncertainties are defined by 
constructing an envelope that embraces the three separate PDF sets from the respective groups, 
together with their individual associated uncertainties. Within each group, PDF reweighting 
1^0 1 is used to evaluate the respective uncertainties. Additional uncertainties from the depen- 
dence on the strong coupling constant are estimated with MSTW08 PDF sets. The resulting 
uncertainty on the integrated SM DY distribution for masses above 1.1 TeV, from all uncertain- 
ties related to the choice of PDF, is estimated to be 13%. 



CMS \/i=7TeV, £ =2.3fb ' CMS \/i"= 7 TeV, /3= 2.1 fb ' 




M,,[GeV] M,-e.[GeV] 

Figure 1: Dimuon (left) and dielectron (right) invariant mass spectra compared with the SM 
predictions and a simulated ADD signal with Aj = 2.8 TeV (ADD K-factor 1.0, no signal trun- 
cation). The highest-mass bins contain all contributions above 2.3 TeV. The error bars reflect 
the statistical uncertainty. 

Contributions from tt, tW, and EW vector boson pair production to the dimuon and dielectron 
mass spectrum are estimated by using simulations with MadGraph [41] and PYTHIA 6. The 
background contributions are cross-checked with a control sample dominated by these pro- 
cesses, including events with an electron and a muon passing requirements similar to those 
used for the signal leptons. Taking into account the differences in the acceptance and efficien- 
cies between muons and electrons, the ratios between the expected ee, fiji, and efi backgrounds 
from the tt, tW, and diboson contributions in the SM are well understood from lepton univer- 
sality. The measured efi mass spectrum is found to be well reproduced by the simulations. 
The agreement has been confirmed up to masses of ~ 500 GeV, above which the statistical 
uncertainties on the e}i spectrum become large. 

Background contributions at large dimuon mass from multijet processes and cosmic-ray muons 
are negligible for our event selection requirements. 

In addition to those backgrounds that are common with the dimuon channel, the dielectron 
channel receives background contributions from multijet events with 2 jets that pass the elec- 
tron selection and W + jets events with 1 jet passing electron selection. Events of the type 
7 + jets, where the photon converts to e+e^ and both the photon and a jet are reconstructed as 
electrons that pass selection, are also considered. The rate for jets to be reconstructed as elec- 
trons is determined from a control sample of events selected by a single-electromagnetic-cluster 
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trigger with a lower threshold. The electron selection criteria, including the isolation require- 
ments, are relaxed to define electron candidates in this sample. Events are required to have 
no more than 1 such reconstructed electron to suppress the contribution from the DY process. 
Residual contributions from W + jets and 7 + jets events in the control sample are estimated 
from simulation. The estimated probability for an electron candidate to pass the full set of elec- 
tron selection criteria is then used to weight the events that have 2 such candidates passing the 
double-electromagnetic-cluster trigger. 

Both for dimuon and dielectron events, the contributions from non-DY processes sum up to 
less than 10% of the expected backgound in the signal region. 

Using Z-candidate events, detailed studies are performed of possible differences in the elec- 
tron and muon reconstruction efficiencies of simulated events and data [42]. No statistically 
significant deviations between data and simulations are found, indicating that the simulated 
lepton reconstruction efficiencies are reliable. In both channels, uncertainties on the simulated 
acceptance and reconstruction efficiencies at large dilepton mass are included in the statisti- 
cal evaluation of the result. Uncertainties related to momentum reconstruction of muons and 
energy estimation of electrons are also taken into account. 

The systematic uncertainties for the integrated dimuon and dielectron backgrounds in the sig- 
nal region are summarized in Table [l] With the exception of the uncertainty on the integrated 
luminosity, which is treated as fully correlated between the 2 channels, they are assumed to be 
independent. 

Figure [l] shows the observed and expected mass distributions in the 2 search channels as a 
function of dilepton mass. Measurements and SM predictions are found to be in agreement 
within statistical and systematic uncertainties. In both channels, no significant excess of events 
is observed in the high-mass region, and no events are found in the signal region. The cor- 
responding SM background expectation for dilepton mass above Mmin = 1-1 TeV is 1.0 ± 0.2 
events in the dielectron channel and 1.3 ± 0.2 events in the dimuon channel. The observed 
number of events N^bs and the SM expectation are in agreement in several control regions, as 
shown in in Table IH 

Table 1: Summary of systematic uncertainties for the integrated dimuon and dielectron invari- 
ant mass spectra in the signal regions. 



Systematic uncertainty 


Uncertainty 
on signal (%) 


Uncertainty 
on background (%) 


Integrated luminosity 

Trigger and reconstruction efficiency 

Muon momentum resolution 

Electron energy scale 

Drell-Yan PDF uncertainties 

Drell-Yan higher order corrections 


4.5 

1 

1-3 


4.5 

3 {}l^),3{ee) 
5 

1-3 
13 
10 



For the statistical evaluation of the measurements, we count the observed events above Mmin- 
For each channel, the probability of observing Nobs events in the signal region is given by a 
Poisson distribution. The statistical model for the Poisson means includes parameters that are 
used to describe the influence of the systematic uncertainties listed in Table [l] on the expected 
signal and background events. Limits on the cross sections for signals in the regions of accep- 
tance are calculated with a CLs approach ||43| . The applied test statistic is a one-sided profile 
likelihood ratio [44J corresponding to the selected models. The systematic uncertainties are 
included in the statistical evaluation by extending the likelihood with additional probability 



6 



Table 2: Comparison of the observed and expected number of events in control and signal 
regions for the dimuon and dielectron mass distributions. Expected signal contributions are 
shown for Aj = 2.8 TeV (ADD K-factor 1.0, signal truncation at Mmax = At). 



Hpi, C = 2.3 fb"^ 


ee, £ = 2.1 fb"^ 


Mass 


Nobs 


Background 


Signal exp. 


Mass 


Nobs 


Background 


Signal exp. 


region [TeV] 




expectation 


At = 2.8 TeV 


region [TeV] 




expectation 


At = 2.8 TeV 


Control regions 


Control regions 


0.14-0.20 


3723 


3690±300 




0.12-0.20 


6592 


6598 ±530 




0.20-0.40 


1674 


1605±160 




0.20-0.40 


1413 


1301±120 




0.40-0.60 


131 


122±13 




0.40-0.60 


88 


103±11 




0.60-0.80 


16 


21±3 




0.60-0.80 


21 


18±3 




0.80-1.10 


8 


5±1 


0.8 


0.80-1.10 


7 


6±1 


0.6 


Signal region 


Signal region 


> 1.10 





1.0±0.2 


3.2 


> 1.10 





1.3±0.2 


2.7 



density functions that parameterize the respective uncertainties. The exclusion threshold is set 
to CLs = 0.05 ( > 95% confidence). 

The RooStats |45| software for statistical data analysis is used for the numerical evaluation of 
the CLs limits. At 95% confidence level (CL), we exclude signal cross sections above 1.2 fb 
(1.8 fb expected) in the dimuon channel and 1.6 fb (2.3 fb expected) in the dielectron channel. 
The combined upper limit at 95% CL on the signal cross section in both channels crs,ppt+ee is 
found to be 1.4 fb, while the expected limit is 2.2 fb. 

The observed limits on tJg are translated into exclusion limits on the ADD parameters. To ac- 
count for interference effects, the expected signal contribution for a particular choice of model 
parameters is evaluated by subtracting the SM DY cross section at LO from the cross section 
with the ADD LO contributions. Limits are based either on the leading-order ADD scenario 
without higher-order corrections or on an assumed higher-order correction factor (K-factor) of 
1.3 for the ADD signal contributions. Based on studies of QCD NLO corrections to dilepton 
processes in the ADD model Il46ll47l , the K-factor of 1.3 corresponds to a conservative choice. 
Figure |2] shows the limits for the HLZ convention for different ranges of validity of the model, 
assuming a K-factor of 1.3 and no signal contribution beyond the cutoff Mmax- Table |3] summa- 
rizes the limits on the GRW and HLZ parameters for truncation of the signal at Mmax = or 
Mmax = Mg. Results are also given for an evaluation of limits separately for dimuon or dielec- 
tron measurements. Including our recently published results on diphoton events [19], which 
have comparable sensitivity, improves the observed combined limits on Ms presented in Table 
|3|by 0.1 (0.3) TeV for n = 2 and 0.1 (0.1) TeV forn > 3 without (with) K-factors for the signal 
contributions. 

In summary, a search for the effects of large extra dimensions in dimuon and dielectron invari- 
ant mass spectra using the CMS detector at the LHC has been presented. The results are found 
to be in agreement with SM predictions, and no significant excess of events is observed at large 
values of dimuon or dielectron mass. In the region of dilepton masses above 1.1 TeV, no events 
are found. Our results extend the limits on ADD models based on the analysis of dilepton sig- 
natures. The combination with diphoton results provides the most stringent limits on graviton 
decay in the ADD framework to date. 

We thank M.C. Kumar, R Mathews, and V. Ravindran for useful discussions on QCD NLO 
corrections in the ADD model. We wish to congratulate our colleagues in the CERN acceler- 
ator departments for the excellent performance of the LHC machine. We thank the technical 
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Figure 2: Observed and expected 95% CL lower limits on Ms, obtained by combining the fiji 
ee results, for different numbers of extra dimensions n, applying a signal K-factor of 1.3. A 
confidence interval for the expected limit corresponding to 2 standard deviations (sd) is shown 
for the case n = 3. 



Table 3: Observed lower limits in TeV at 95% CL within GRW and HLZ conventions for trun- 
cation at Mmax = At (GRW) or M^ax = Ms (HLZ). 



ADD K-factor 


At [TeV] (GRW) 


Ms [TeV] (HLZ) 
n = 2 n = 3 n=4 n = 5 n = 6 n = 7 


HH, o-s,mi< 1.2 fb (1.8 fb expected) at 95% CL 


1.0 


2.8 


3.0 3.4 2.8 2.5 2.3 2.2 


1.3 


3.0 


3.2 3.5 3.0 2.7 2.4 2.3 


ee, (7s,ee< 1-6 fb (2.3 fb expected) at 95% CL 


1.0 


2.8 


2.9 3.3 2.8 2.5 2.3 2.2 


1.3 


2.9 


3.1 3.4 2.9 2.5 2.4 2.2 


piji and ee, crs„^^+ee< 1-4 fb (2.2 fb expected) at 95% CL 


1.0 


3.1 


3.7 3.7 3.1 2.8 2.5 2.4 


1.3 


3.2 


3.8 3.8 3.2 2.9 2.7 2.5 
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